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ly^ ■ Abstract: The apphcation of /cf-factorization supplemented with the CCFM small- 

. X evolution equation to heavy quark production at the TEVATRON and at HERA is 

^ . discussed. The bb production cross sections at the TEVATRON can be consistently 

\ described using the /cj-factorization formalism together with the unintegrated gluon 

^ \ density obtained within the CCFM evolution approach from a fit to HERA F2 data. 

O ' Special attention is drawn to the comparison with measured visible cross sections, 

^ ■ which are compared to the hadron level Monte Carlo generator Cascade. 



^ ■ 1 Introduction 

The calculation of inclusive quantities, like the structure function F2{x,Q'^) at HERA, per- 
1-^ ' formed in NLO QCD is in perfect agreement with the measurements. However, Catani argues, 
^ ■ that the NLO approach, although phenomenologically successful for F2{x, Q"^), is not fully satis- 
^ . factory from a theoretical viewpoint, because "the truncation of the splitting functions at a fixed 
perturbative order is equivalent to assuming that the dominant dynamical mechanism leading 
to scaling violations is the evolution of parton cascades with strongly ordered transverse mo- 
menta" [ As soon as exclusive quantities like jet or heavy quark production are investigated, 
the agreement between NLO coefficient functions convoluted with NLO DGLAP [ H, ^ |^ 
parton densities and the data is not at all satisfactory: large so-called i^'-factors (normalization 
factors) [ 1^, 0, H, 1^ are needed to bring the NLO calculations close to the data {K ~ 2 — 4 for 
bottom production at the TEVATRON), indicating that in the calculations a significant part 
of the cross section is still missing. 

At small X the structure function F2{x, Q"^) is proportional to the sea quark density, and the 
sea-quarks are driven via the DGLAP evolution equations by the gluon density. Standard QCD 
fits determine the parameters of the initial parton distributions at a starting scale Qq. With 
help of the DGLAP evolution equations these parton distributions are then evolved to any other 
scale Q^, with the splitting functions still truncated at fixed 0(a<j) (LO) or 0{a1) (NLO). Any 
physics process in the fixed order scheme is then calculated via coUinear factorization into the 
coefficient functions C"(-) and coUinear (independent of kt) parton density functions: fa{z, Q"^): 
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At large energies (small x) the evolution of parton densities proceeds over a large region in 
rapidity Ay ~ log(l/a;) and effects of finite transverse momenta of the partons may become 
increasingly important. Cross sections can then be kt - factorized [ ^ into an off-shell {kt 
dependent) partonic cross section a{^,kt) and a kt - unintegrated parton density function 
T{z, kt): 
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The unintegrated gluon density J-'{z, kt) is described by the BFKL [ ^ |T2|, |T3[ evolution 
equation in the region of asymptotically large energies (small x). An appropriate description 
valid for both small and large x is given by the CCFM evolution equation [ |1^, |1^, [1^, |T^, 



resulting in an unintegrated gluon density A{x, kt, g), which is a function also of the additional 
evolution scale q described below. 

In [ |18], [l|] Catani argues that by explicitly carrying out the kt integration in eq.(|^) one can 



obtain a form fully consistent with collinear factorization: the coefficient functions and also the 
DGLAP splitting functions leading to fa{z,Q'^) are no longer evaluated in fixed order pertur- 
bation theory but supplemented with the all-order resummation of the as log 1/x contribution 
at small x. 

The application of kt - factorization to heavy quark hadroproduction is discussed in detail 
in [ ig, ^ In the present paper bottom production at the TEVATRON and at HERA 
is investigated using the kt - factorization approach. The unintegrated gluon density has been 
obtained previously in from a CCFM fit to the HERA structure function ^2(0;, Q^). All 
free parameters are thus fixed and absolute predictions for bottom production can be made. 
The universality of the unintegrated CCFM gluon distribution will be tested by the application 
to TEVATRON and HERA results. 

The basic features of the CCFM evolution equation are recalled and the unintegrated gluon 
density is investigated. Then the calculations for bb production at the TEVATRON is presented 
as well as calculations of the visible cross section for hb production at HERA. 



2 The CCFM evolution equation 

A solution of the CCFM evolution equation, which properly describes the inclusive structure 
function F2{x, Q^) and also typical small x final state processes at HERA has been presented in 
detail in [ Figure |l| illustrates the pattern of QCD initial-state radiation in a small-x event 
at a pp collider, together with labels for the kinematics. According to the CCFM evolution 
equation, the emission of partons during the initial cascade is only allowed in an angular-ordered 
region of phase space. The maximum allowed angle S is defined by the hard scattering quark 
box, producing the heavy quark pair. In terms of Sudakov variables the quark pair momentum 
is written as: 

p,+p,-=T(p« + Sp(2)) + Qi (3) 

where p^^^ and p'-^-' are the four- vectors of incoming protons, respectively and Qt is the transverse 
momentum of the quark pair in the laboratory frame. Similarly, the momenta pi of the gluons 
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Figure 1: Kinematic variables for multi-gluon emission. The t-channel gluon four-vectors are 
given by ki and the gluons emitted in the initial state cascade have four-vectors Pi. The minimum 
(maximum) angle for any emission is obtained from the quark box, as indicated with S. 



emitted during the initial state cascade are given by (here treated massless): 



= V^{P^'^ + ^^P^'^) + PU , ir=^. (4) 

svt 

with Vi = (1 — Zi)xi-i, Xi = ZiXi-i and s = (p^^^ being the squared center of mass energy. 

The variable C,i is connected to the angle of the emitted gluon with respect to the incoming 
proton and Xi and Vi are the momentum fractions of the exchanged and emitted gluons, while 
Zi is the momentum fraction in the branching (i — 1) — >■ i and pu is the transverse momentum 
of the emitted gluon i. 

The angular-ordered region is then specified by (for the lower part of the cascade in Fig. 
the upper part is obtained by properly exchanging the variables): 

^0 < 6 < ■ • • < < S (5) 

which becomes: 

Zi-iQi-i < Qi (6) 
where the rescaled transverse momenta qi of the emitted gluons is defined by: 

qi = ^/s^i = (7) 
i Zi 

The CCFM equation for the unintegrated gluon density can be written [ |l^ ^ ^ ^ as 
an integral equation: 

/dz f d?'q ~ /x\ 

~ — - ^(l)^s{q,zq)P{z,q,kt)A{-,k[,q] (8) 
z J T^q ^ z ' 
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with k[ = \kt + {l — z)q\ and q being the upper scale for the last angle of the emission: q > Znqn, 
qn > Zn-iqn-i, li > Qo- Here q is used as a shorthand notation for the 2-dimensional vector 
of the rescaled transverse momentum q = qt= p^/ [1 — z). The splitting function P{z, q, kt) and 
the Sudakov form factor As{q, zq) are given explicitly in [ |22| . 



2.1 The unintegrated gluon density 



In [|2^ the unintegrated gluon density xA{x, kf, q) has been obtained from a fit to the structure 
function F2{x,Q'^) Q 

In Fig. 1^ the CCFM unintegrated gluon density distribution as a function of x and k'f is 
shown and compared to 

dxG{x, jj? 



J-{x, k^ 
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with xG{x,^) being the collinear gluon density of GRV 98 [ |2B| in LO and NLO. 



The unintegrated gluon density can be related to the integrated one by: 

^G{x,fi)\ / dklxA{x,kl,q) (10) 

Here the dependence on the scale of the maximum angle q is made explicit: the evolution 
proceeds up to a maximum angle related to g, which plays the role of the evolution scale in the 
collinear parton densities. This becomes obvious since 

f = xfEs = x^^^xfh = s + Ql (11) 

The last expression is derived by using Pq + pq — xf^p2 + x^ppi + Qt-, S ~ x^P /xf^ and 

and s = x^x^ s — Qf. This can be compared to a possible choice of the renormalization and 
factorization scale /i^ in the collinear approach with /i^ = + 4 ■ itiq and the similarity between 
H and q becomes obvious. 

In Fig. 1^ the CCFM gluon density integrated over kt according to eq.(|l^) is compared to the 
gluon densities of GRV 98 [ ^ in LO and NLO. It is interesting to note that the CCFM gluon 
density is flat for a; at the input scale Q = 1 GeV. Even at larger scales the collinear gluon 
densities rise faster with decreasing x than the CCFM gluon density. However, after evolution 
and convolution with the off-shell matrix element the scaling violations of F2(x,Q^) and the 



rise of F2 towards small x is reproduced, as shown in [ |2^, Fig. 4 therein]. A similar trend is 
observed in the collinear fixed order calculations, when going from LO to NLO: at NLO the 
gluon density is less steep at small x, because part of the x dependence is already included in 
the NLO Pqg splitting function, as argued in [|l|. 

From Fig. |^ one can see, that the integrated gluon density from CCFM is larger in the 
medium x range, than the ones from the collinear approach. Due to an additional l/k'j sup- 
pression in the off-shell matrix elements, the gluon density obviously has to be larger to still 
reproduce the same cross section. In addition only gluon ladders are considered in the k^ - 
factorization approach used here, which means that the sea quark contribution to the structure 
function F2{x,Q'^) comes entirely from boson gluon fusion, without any contribution from the 



intrinsic quark sea. One also should remember, that the relation in eq.(lO) is only approxi- 
mately true, since the gluon density itself is not a physical observable. 

Fortran program for the unintegrated gluon density xA{x, , q) can be obtained from [ Eq] 
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Figure 2: The CCFM kf dependent (unintegrated) gluon density [^^, |^/ at q = 10 GeV as a 
function of x for different values of ( upper) and as a function of k^ for different values of 
X (lower) compared to '^^'^^^2'^ ^ with xG{x,n) being the collinear gluon density of GRV 98 [\2dij 
in LO and NLO. 



3 bb production at the TEVATRON 



The cross section for bb production in pp collision at ^/s = 1800 GeV is calculated with Cas- 
[ |2^, which is a Monte Carlo implementation of the CCFM approach described 
The off-shell matrix element as given in [ Il0| for heavy quarks is used with mi, = 4.75 



CADE 

above 



GeV. The scale /i used in ^^(yU ) is set to fi 



mi + p^ (as in 



with Pt being the 



transverse momentum of the heavy quarks in the pp center-of-mass frame. In Fig. |^ the predic- 
tion for the cross section for bb production with pseudo-rapidity \y^\ < 1 is shown as a function 
of p™" and compared to the measurement of DO [ . Also shown is the NLO prediction from [ 
p7| , taken from [ ^]. In Fig. ^ the measured cross section of CDF [ ^ is shown for 3 values of 
with the kinematic constraint of \y'^\^ \y^\ < 1 and p^^"{b) > 6.5 GeV together with the 
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Figure 3: The CCFM gluon density (solid line) integrated over kt as a function of x for different 
values of Q. For comparison the CRV 98 [ gluon density in LO (dashed line) and NLO 
(dotted line) is also shown. 



prediction from CASCADE and the NLO calculation from [ taken from [ ^]. In all cases the 
NLO calculation used mf, = 4.75 GeV and the factorization and renormalization scales were set 
+ p%. Both, DO and CDF measurements are above the NLO predictions by a 



m] 



to ji = ruq 

factor of ~ 2. The Cascade predictions are in reasonable agreement with the measurements. 
A similarly good description of the DO and CDF data has been obtained in [ |3^ using also kt - 
factorization but supplemented with a BFKL type unintegrated gluon density. It is interesting 
to note, that the CCFM unintegrated gluon density has been obtained from inclusive F2{x, Q^) 
at HERA. In this sense, the prediction of Cascade is a parameter free prediction of the bb 
cross section in pp collisions. This also shows for the first time evidence for the universality of 
the unintegrated CCFM gluon distribution. 

Since Cascade generates full hadron level events, direct comparisons with measured cross 
sections for the production of b quarks decaying semi-leptonically into muons are possible. The 
muon cross sections as a function of the transverse momentum pi^ and pseudo-rapidity \y^\ as 



measured by DO [ ^ are compared to the Cascade prediction in Fig. ^ and Fig. |^. Both 
the and {y'^l cross sections are well described, whereas NLO calculations underestimate the 
cross sections by a factor of ~ 4 in the small angle range [ 2^ . 

It has been argued in [ that within the coUinear approach supplemented with DGLAP 
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Figure 4: Cross section for bb production 
with \y^\ < 1 as a function of p^^"^^ . Shown 
are the DO [ 0/ data points, the fixed or- 
der NLO prediction, and the prediction of 
Cascade. 



■ CDF Data, L = 17.4pb'' 
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mi=4.75 GeV/c^ A5 = 1 40 MeV 
/i=Mo=((PT/+P,/)/2 + m/)"' 
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Figure 5: Cross section for bb production 
with \y^\ < 1 and p'!^^"'{b) > 6.5 CeV as a 
function of p^''"'(b). Shown are the CDF [ 
^ data points, the fixed order NLO predic- 
tion, and the prediction o/ CASCADE. 
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Figure 6: Cross section for muons from b- 
quark decays as a function of pi^ (per unit 
rapidity) as measured by DO compared 
to the prediction of CASCADE . 
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Figure 7: Cross section for muons from b- 
quark decays as a function of |?/^| for two 
different cuts as measured by DO [ p^/ 
compared to the prediction of Cascade . 



parton showers, heavy quark excitation in form of the processes Qg —>■ Qg and Qq Qq con- 
tributes significantly to the bb cross section in the central rapidity region. The bb measurements 
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at the TEVATRON can be reasonably well described [ ^| by the price of including heavy quark 
excitation and a heavy quark component in the structure function of the proton. However kt- 
factorization is superior to the collinear approach with heavy quark excitation, since in the 
/crfactorization approach heavy quarks are produced perturbatively only via hard scattering 
matrix elements, without further additional assumptions. 



X 




> 

(U 

o 



0.1 



0,08 



X3 



0.06 - 



.0.04 



0.02 





CASCADE 




P^THIA 







1 



[GeV'] 



Figure 8: Comparison of Xn and kfn distributions of the gluons entering the process g+g b+b. 
Shown are the predictions from CASCADE representing kffactorization with the CCFM unin- 
tegrated gluon density and also from Pythia representing the collinear approach supplemented 
with initial and final state DGLAP parton showers to cover the phase space of pt ordered QCD 
cascades. 



In Fig. ^ the x„ and ktn distributions of the gluons entering the hard scattering process 
g + g ^ b -\-b (see Fig. |I]) are shown and the predictions from the kt - factorization approach 
(Cascade ) are compared to the standard collinear approach (here Pythia [ ^ with LO 
gg bb matrix elements supplemented with DGLAP parton showers to simulate higher order 
effects). Whereas the x„ distributions agree reasonably well, a significant difference is observed 
in the ktn distribution. However this is not surprising: the kf factorization approach includes 
a large part of the fixed order NLO corrections (in collinear factorization). Such corrections 
are gg — >■ QQg, where the final state gluon can have any kinematically allowed transverse 
momentum, which could be regarded as a first step toward a non-p^r ordered QCD cascade. 
The importance of higher order QCD radiation can be estimated from the number of hard gluon 
radiation Ng with a transverse momentum of the gluon > m^, with mx = (Pr + Pr)/2 + ^b, 
by defining the probability Ph- 



(12) 



In a leading order calculation in the collinear approach, even supplemented with DGLAP 
parton showers, it is expected that Ph{Ng > 1) — > 0, because of the approximate transverse 
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Figure 9: Probability Ph{Ng 
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of hard gluon radiation with pf > mx in bb production. 



Shown are calculations in the kt - factorization approach (^Cascade^ and from the LO DGLAP 
approach fPYTHlAj. 



ordering constraint in DGLAP. In a fixed order a1 (NLO) calculation, Ph{Ng = 1) 7^ 0, but 
Ph{Ng > 2) = 0. In Fig. |^ the probability Ph is shown as a function of Ng for the production of 
bb events with p^ > 5 GeV, \y''\ < 1 obtained in the kt - factorization approach via Cascade. 
Also shown for comparison is the collinear approach in LO supplemented with DGLAP parton 
showers obtained from Pythia. In the fc^-factorization approach a significant contribution of 
hard higher order QCD radiation with p^ > rriT is obtained. 



4 bb production at HERA 



In [ ^ the prediction of Cascade for the total bb cross section was compared to the extrap- 
olated measurements of the HI [ ^ and ZEUS [ Q experiments at HERA. Since Cascade 
generates full hadron level events, a direct comparison with measurements can be done, be- 
fore extrapolating the measurement over the full phase space to the total bb cross section. 
ZEUS [ 1^ has measured the dijet cross section which can be attributed to bottom production 
by demanding an electron inside one of the jets. In the kinematic range of < 1 GeV^, 
0.2 < y < 0.8, at least two jets with E^^^^"^^ > 7(6) GeV and \ri^'^^\ < 2.4 and a prompt electron 
with p^ > 1.6 GeV and \ri^ \ < 1.1, ZEUS [|| quotes the cross section: 

a,V,!:e++d.,e*+e-+x = 24.9 ± QAt^i pb (ZEUS [ i) (13) 

with the statistical (first) and systematic (second) error given. Within the same kinematic 
region and applying the same jet algorithm Cascade predicts: 

a,V;_:e++d»,et+e-+x = 20.3^1;^ pb, (14) 
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where the error reflects the variation of nif, 



measurement within the statistical error. It has been shown in [ ^ 



4.75 =F 0.25 GeV. This value agrees with the 

that the extrapolation 



to the bb cross section cross section includes large model uncertainties and the comparison of 
extrapolated cross sections with predictions from NLO calculations are questionable. 
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Figure 10: The differential cross section da /dx?^^ for heavy quark decays as measured by ZEUS [ 
^ and compared to CASCADE. Shown are the charm, bottom and the sum of both contributions 
to the cross section (Note there is no additional K factor applied). 
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It is interesting to note, that the Cascade result agrees well with the Pythia result for 
the di-jet plus electron cross section, if heavy quark excitation is included in Pythia. As in 
the case of bb production at the TEVATRON higher order QCD effects are important, which 
are already included in the kt - factorization approach. 

In Fig. [ly the differential cross section for heavy quark decays (charm and bottom) as a 
function of predicted by Cascade is compared with the measurement of ZEUS [0. A 
significant fraction of the cross section has < 1, which is similar to the observation in 



charm photo-production [22]. In LO in the collinear factorization approach this is attributed 
to resolved photon processes. However, in kf - factorization, the x^ distribution is explained 
naturally [ ^ , because gluons in the initial state need not to be radiated in a pt ordered region 



and therefore can give rise to a high px jet with transverse momentum larger than that of the 
heavy quarks. 

The prediction of Cascade has also been compared to the measurement of HI [ |^ for 
electro-production cross section in < 1 GeV^, 0.1 < y < 0.8, < 2 GeV and 35° < 6"^ < 
130°: 

a{ep e'bbX fiX') = 0.176 ± 0.016{stat.)tom7isyst.) nb 
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This cross section already includes the extrapolation from measured jets to the muon. In the 
same kinematic range CASCADE predicts: 



a{ep e'bbX ^X') 



0.0661°:°°? nb. 



which is a factor ~ 2.6 below the measurement. It has been shown in that the experimental 
results from HI and ZEUS are different by a factor ~ 2. It has been checked, that the difference 
between the experiments is not due to different kinematics. 



5 Conclusion 

Bottom production at the TEVATRON can be reasonably well described using the kt- factor- 
ization approach with off-shell matrix elements for the hard scattering process. Also the cross 
section for b quarks decaying semi-leptonically into muons at small angles is reasonably well 
described, whereas NLO calculations fall below the data by a factor of ~ 4. One essential in- 
gredient for the satisfactory description is the unintegrated gluon density, which was obtained 
by a CCFM evolution fitted to structure function data at HERA, showing evidence for the 
universality of the unintegrated gluon density. The comparison with the data was performed 
with the Cascade Monte Carlo event generator, which implements /cj-factorization together 
with the CCFM unintegrated gluon density. 

Measurements of bottom production at HERA are also compared to predictions from Cas- 
cade. The visible dijet plus electron cross section attributed to 6-production as measured 
by ZEUS could be reproduced within the statistical error. However, the situation is different 
with the HI measurement: the visible muon cross section is already a factor of 2.6 above the 
prediction from CASCADE. Further measurements, also differential, are desirable to clarify the 
situation at HERA. 

In general the fcf-factorization approach has now proven to be successful in a wide kinematic 
range. It is worthwhile to note, that the approach presented here, in addition to the description 
of bb production, is also able to describe charm and jet-production at HERA as well as the 
inclusive structure function F2. It is the advantage of this approach that contributions to the 
cross section, which are of NLO or even NNLO nature in the coUinear ansatz, are consistently 
included in fcf-factorization due to the off-shellness of the gluons, which enter into the hard 
scattering process. 



6 Acknowledgments 

I have learned a lot from discussions with B. Andersson, G. Gustafson, G. Ingelman, L. Lonnblad 
and T. Sjostrand. I am grateful to G. Salam for all his patience and his help in all different 
kinds of discussions concerning CCFM and a backward evolution approach. Many thanks also 
go to S. Baranov and N. Zotov who were the first raising my interest in fc^-factorization. I am 
grateful to E. Elsen, L. Gladilin, G. Ingelman, L. Jonsson, P. Schleper for careful reading and 
comments on the manuscript. No words are adequate for everything I have with Antje. I also 
want to thank the DESY directorate for hospitality and support. 



11 



References 



[1] S. Catani, Aspects of QCD, from the Tevatron to LHC, in Proceedings of the Inter- 
national Workshop Physics at TeV Colhders (Les Houches, France, 8-18 June, 1999), 
hep-ph/0005233. 

[2] V. Gribov, L. Lipatov, Sov. J. Nucl. Phys. 15(1972)438 and 675. 

[3] L. Lipatov, Sov. J. Nucl. Phys. 20(1975)94. 

[4] G. Altarelh, G. Parisi, Nucl. Phys. B 126(1977)298. 

[5] Y. Dokshitser, Sov. Phys. JETP 46 (1977) 641. 

[6] CDF Collaboration; F. Abe et al., Phys. Rev. D 55(1997)2546. 

[7] DO Collaboration; B. Abbott et al., Phys. Lett. B 487(2000)264. 

[8] HI Collaboration; C. Adloff et al., Phys. Lett. B 467(1999) 156, and erratum ibid. 

[9] ZEUS Collaboration; J. Breitweg et al, Eur. Phys. J. C 18(2001)625. 
[10] S. Catani, M. Ciafaloni, F. Hautmann, Nucl. Phys. B 366(1991) 135. 
[11] E. Kuraev, L. Lipatov, V. Fadin, Sov. Phys. JETP 44 (1976) 443. 
[12] E. Kuraev, L. Lipatov, V. Fadin, Sov. Phys. JETP 45 (1977) 199. 
[13] Y. Balitskii, L. Lipatov, Sov. J. Nucl. Phys. 28(1978)822. 
[14] M. Ciafaloni, Nucl. Phys. B 296 (1988) 49. 

[15] S. Catani, F. Fiorani, G. Marchesini, Phys. Lett. B 234(1990)339. 
[16] S. Catani, F. Fiorani, G. Marchesini, Nucl. Phys. B 336 (1990) 18. 
[17] G. Marchesini, Nucl. Phys. B 445(1995)49. 

[18] S. Catani, fcf-factorisation and perturbative invariants at small x, in Proceedings of the 
International Workshop on Deep Inelastic Scattering, DIS 96 (Rome, Italy, 15-19 April, 
1996), hep-ph/9608310. 

[19] E. M. Levin, M. G. Ryskin, Y. M. Shabelski, A. G. Shuvaev, Sov. J. Nucl. Phys. 54(1991) 
867. 

[20] M. G. Ryskin, A. G. Shuvaev, Y. M. Shabelski, Phys. Atom. Nucl. 64(2001) 120, 
|arXiv:hep-ph/9907507 . 

[21] M. G. Ryskin, A. G. Shuvaev, Y. M. Shabelski, Comparison of /cT^-factorization approach 
and QCD parton model for charm and beauty hadroproduction, 2000, [hep-ph/0007238| . 

[22] H. Jung, G. Salam, Eur. Phys. J. C 19 (2001) 351, hep-ph/0012143. 

[23 



G. Bottazzi, G. Marchesini, G. Salam, M. Scorletti, JHEP 12 ( 1998) Oil, |h^-ph/981054e . 

12 



[24] J. Kwiecinski, A. Martin, P. Sutton, Phys. Rev. D 52 (1995) 1445. 

[25] H. Jung, The CCFM Monte Carlo generator Cascade for lepton - proton and proton - pro- 
ton collisions, Lund University, accepted by Comp. Phys. Comm., 2001, hep-ph/0109102, 
DESY 01-114, http://www.quark.lu.se/~hannes/cascade/. 

[26] M. Gliick, E. Reya, A. Vogt, Eur. Phys. J. C 5 (1998) 461, hep-ph/9806404. 

[27] S. Frixione, M. Mangano, Nucl. Phys. B 483(1997)321. 

[28] M. Mangano, P. Nason, G. Ridolfi, Nucl. Phys. B 373 (1992) 295. 

[29] DO Collaboration; B. Abbott ct al., Phys. Rev. Lett. 84(2000)5478. 

[30] P. Hagler et al., Phys. Rev. D62 (2000) 071502. 

[31] E. Norrbinn, T. Sjostrand, Eur. Phys. J. C 17(2000)q37, hep-ph/0005110. 
[32] R. Field, private communication (unpublished). 
[33] T. Sjostrand, Comp. Phys. Comm. 82 (1994) 74. 

[34] H. Jung, Unintegrated parton densities appHed to heavy quark production in the CCFM 

approach, in Proceedings of the Rinherg workshop on "New trends in HERA physics", 
Ringberg Castle, Tegernsee, Germany. (2001), hep-ph/0109146. 

[35] S. Baranov, N. Zotov, Phys. Lett. B 491(2000) 111. 



13 



